Remote Dipolar Interactions for Objective Density Calibration and Flow Control of 

Excitonic Fluids 
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In this paper we suggest a method to observe remote interactions of spatially separated dipolar 
quantum fluids, and in particular of dipolar excitons in GaAs bilayer based devices. The method 
utilizes the static electric dipole moment of trapped dipolar fluids to induce a local potential change 
on spatially separated test dipoles. We show that such an interaction can be used for a model- 
independent, objective fluid density measurements, an outstanding problem in this field of research, 
as well as for inter-fluid exciton flow control and trapping. For a demonstration of the effects on 
realistic devices, we use a full two-dimensional hydrodynamical model. 



Dipolar quantum fluids are very intriguing physical en- 
tities that are intensively sought for in various material 
systems, such as ultra-cold polar molecules h| and ex- 
citons in semiconductor bilayer structures UQ. The 
interest in such fluids arises from the basic unresolved 
questions regarding their possible rich thermodynamic 
phases, which are predicted to emerge from the inter- 
play between their quantum statistics and the particular 
form of the dipole-dipole interaction 0-0 • While the 
local interactions between nearby dipoles is expected to 
determine the physical state of the fluid, it should be 
interesting to observe the effect of interactions between 
spatially remote fluids. This concept of remote inter- 
actions can also be utilized, as we will suggest here, in 
a couple of aspects: it can be used to determine in an 
objective, model-independent way the density of a dipo- 
lar fluid, a tricky "boot-strap" problem, especially for 
excitons in bilayer systems. It can also be utilized as a 
method to control the flow of one fluid by other, remotely 
located fluids. 

As mentioned above, dipolar fluids can be realized for 
excitons in bilayer structures. An exciton consists of a 
pair of an electron and a hole inside a semiconductor, 
which are bound together by their mutual coulomb at- 
traction, forming a boson-like quasi particle. In elec- 
trically gated GaAs coupled quantum wells (CQW), also 
known as bilayer structures, an externally applied electric 
field, using semi-transparent electrical gates, results in 
the creation of two-dimensional dipolar excitons. These 
light mass excitonic quasi-particles have their electron 
and hole constituents separated into two layers (Fig. 
UK a)) and thus have an increased life time, an elevated 
quantum degeneracy temperature (compared to atoms 
or molecules for example), and all carry a dipole moment 
oriented perpendicular to the QW plane. The two dimen- 
sional dipolar excitons (Xd) interact with each other via 
the repulsive dipole-dipole forces, and therefore exciton 
based systems are model systems for exploring a low di- 
mensional interacting fluid hydrodynamics at quantum 
degeneracy conditions. The Xd can also be transformed 
back into light signals in an almost perfect on/off con- 



trolled fashion by switching the external gates electric 
fields, thus can potentially be used for unique functional 
optoelectronic devices. Various methods can also be used 
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FIG. 1. (a) A Schematic band diagram for a GaAs/AlGaAs 
CQW structure without (right) and with (left) external bias, 
showing the formation of indirect (dipolar) excitons. (b) A 
schematic sketch of the proposed sample structure, contain- 
ing a stack of CQWs (their number is marked by Ncqw), 
a homogeneous back electrode and top semitransparent elec- 
tric gate. The specific top gate geometries for the various 
implementations discussed in the text are given in Fig. [2] 

to effectively transport and manipulate Xd flows. Indeed, 
in the past few years several groups introduced relatively 
efficient ways of controlling Xd fluxes, by using electric 
field gradients Q, surface acoustic waves [H and modu- 
lated gate voltages [9l4l lj . Transport of excitons has also 



been measured in other bilayer systems [12l-ll4j. How- 
ever, all of the above methods are based on external ma- 
nipulations and not on the internal state of the device. 
This may limit future functionality of devices based on 
exciton fluids 0, [ll[ • 

Another crucial experimental difficulty, which is im- 
portant in the aspect of the fundamental understand- 
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ing of the X d fluid thermodynamics, is to independently 
and accurately estimate the density of Xd created in- 
side a sample. Such an estimate is vital for correctly 
identifying and gaining physical insight into phenomena 
such as phase transitions, quantum-degeneracy and par- 
ticle correlations. This task however, turns out to be 
not straight-forward. The current method for estimat- 
ing Xd densities is to deduce it from measurements of 
the energy shift ("blue shift") of recombining X d s pho- 
toluminescence (PL). This energy shift strongly depends 
on the local density correlations around the recombining 
excitons, which in turn depends on the thermodynamic 
state of the fluid 0, [Hj . Therefore this method relies on 
a knowledge of the fluid state, which makes this density 
calibration method a strongly model-dependant, "boot- 
strap" process. Finally, we note that while dipole-dipole 
interactions have been observed to affect the properties 
and the dynamics of a dipolar fluid locally, the remote 
effect of a collection of such dipoles on a system spa- 
tially separated from it has not been observed as far as 
we know. This concept alone is worthwhile pursuing. 

Here, we propose that remote interactions between 
spatially separated excitons can be observed in both spec- 
troscopy and hydrodynamics experiments. Such an ob- 
servation can also be used for resolving practical cur- 
rent challenges as will be suggested in this paper: in the 
first part of the paper we describe a scheme that utilizes 
remote interactions for an objective, model-independent 
dipolar fluid density measurements which does not de- 
pend on local density correlations, and in the second part 
we describe a method to manipulate dipolar fluids using 
their remote interactions. We note that while here we 
present the implementation of the concept for dipolar ex- 
citon fluids in electrostatic traps, the concept developed 
here is quite general and can be implemented in other 
types of trapping devices and also for different types of 
dipolar fluids. 

As a starting point for our discussion, one has to con- 
sider a trap region which will confine the Xd in the QW 
plane, and oppose their spreading due to diffusion and 
repulsive interaction. Various ways for trapping dipo- 
lar excitons have been realized in the past few years 
16|-18|. The approach we took in the current work uti- 
lizes electrostatic traps, which have a flat potential pro- 
file and a sharp boundary Jjl, 2o[, leading to a flat av- 
erage density distribution of excitons in the trap on a 
macroscopic scale, denoted here by n (note that local, 
microscopic density inhomogeneities are always expected 
due to dipole-dipole interactions [6[). These electrostatic 
gates can also define narrow channels, in which excitons 
can propagate only in one direction [21[ . In the following 
discussion, we consider a sample containing a vertical 
stack of separated CQWs (Fig. []Jb)). The number of 
CQW is denoted by Ncqw- As will be shown hereafter, 
increasing Ncqw will increase the effect of the remote 
interactions. We let d denote the dipole length, i.e. the 



spatial separation of the electron and hole wave functions 
in each CQW along the z axis, r is a radius vector in the 
QW plane denoting the position of a test dipole and f x 
is an arbitrary point inside the trapping region. In terms 
of the interaction energy U d (f—f x ) between two dipoles 
separated by a distance \ f — r x \, the total interaction en- 
ergy for a test dipole at a point r is given by |6j : 



E d {f) = N CQ w 



Udif- f x )g d (r- r x )d 2 r x , (1) 



trap 



where g d {r — r x ) is the pair correlation function and 
n9d(r — r x )d 2 r x is the average number of X d s within 
an area d 2 r x at a distance \r — f x \ from the test dipole. 
It is clear that E d (f) for any f inside the trap strongly 
depends on the shape of the correlation function g d . 
Therefore, extracting n from a local measurement of 
E d (f) requires an independent knowledge of g d . This is 
the essence of the objective density calibration problem. 
However, g d (f — r x )\\^_^ x \^ r00 — 1, as at large distances 
all local correlations disappear. Hence, for a test dipole 
remotely located from the rest of the trapped Xd fluid, 
as shown schematically in Fig. Efa), Eq. [T] reduces to: 



E d (f) = N CQ w n U d {r- r x )d 2 r : 



(2) 



trap 



Under such conditions, the dependence of E d {f) on g d 
vanishes and no knowledge of the form of gd is required 
to extract n from a measurement of Ed(f). 

To check this concept of dipolar fluid density mea- 
surements in realistic devices, we consider the geometry 
shown in Fig. [2](b). By applying two different voltages on 
the circular probe gate and the electrostatic ring trap sur- 
rounding it, two potential wells for the dipolar excitons 
in the QW plane are formed, one under each gate. Due 
to the radial gap between the two separated electrodes, 
these potential wells are spatially separated by a narrow 
barrier. Furthermore, as these potential wells have dif- 
ferent depths (due to different applied gate voltages V\ 
and V2), their corresponding Xd PL can be spectrally 
separated [19]. A straight-forward calculation for E d in 
this case using Eq. ([2j) yields the result: 



E d (f =0) = N CQ w -2k 



ne 2 d 2 



Ri R2) ^ 



where R\ and R2 are the inner and outer radii of the 
electrostatic ring trap respectively, and we have used 



U d 



2 d 2 /(e\r- 



), where e is the QW dielectric con- 



stant. This approximation for Ud is valid if \ f — f x \ ^>d 
for every r x , which is the case we are considering here. 



Increasing Ncqw linearly increases the effective de nsity 
at each point f x and thus the interaction strength [221 ]. 
Substituting typical values of Ncqw = 5, n = 10 11 cm -2 , 
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d = 20nm, R\ = 4/ira, R2 = 50/im into Eq. (j3j) gives 
£"d « 0.3raeV. This ^ term can be detected in the PL 
line blue shift of the probe Xd, and can be used to de- 
duce the Xd density n inside the outer trap, as all other 
parameters are known. Since the microscopic density cor- 
relation scale r c is typically smaller or of the order of the 



inter-particle distance ' 



-1/2 



, both Ri and R2 are much 



larger than r c . As a result, this is a model-independent 
density measurement, as it does not assume anything on 
the local state of the Xd fluid in the outer ring. The 
only assumption is that n is constant on the macroscopic 
scale of the trap area, which was already experimentally 



verified [H,[2o|. 
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FIG. 2. (color online), (a) The remote interaction energy 
E at a point r is calculated by integrating the contribution 
from Xd contained inside a specified region, (b) A proposed 
exciton density measurement device. The detected emission 
energy of Xd trapped inside the inner probe gate is shifted 
due to the interaction with Xd in the outer trap. Vi and 
V2 are the applied voltages on the probe and outer trap re- 
spectively, (c) A sketch of a flow device containing a narrow 
channel, surrounded by two half-circular electrostatic traps, 
(d) A numerical calculation of the energy profile along a 2jim- 
wide channel of the device in (c), formed by Xd in the two 
surrounding 10/ira diameter electrostatic traps, which are lo- 
cated lfim away from the channel edges. Here, d = 16nra, 
Ncqw — 4 and n — 10 11 crn~ 2 . 

While in the first part we showed that the effect of 
remote dipolar interactions can be observed spectroscop- 
ically, we now suggest n hydrodynamic effect through 
the concept of flow control using remote interactions. 
We consider a flow device schematically shown in Fig. 
[2fc). Here, a 2/im-wide channel gate is surrounded by 
two half circular electrostatic traps with a 10/im diame- 
ter each. The traps are loaded with a steady state den- 
sity n. Applying Eq. © to all points inside the channel, 
yields an interaction energy profile, shown in Fig. Efd). 
For typical experimental parameters d = 16nm,n = 
10 n cra _2 , Ncqw — 4 and e = ecaAs = 13, the peak has 
a value of « O.lbmeV. We denote this maximal height 
of the energy "bump" by E^. 

We now use an extension of the device shown in Fig. 



[2fc) to demonstrate efficient trapping of Xd in a chan- 
nel. This extension is described in Fig. [3j where the 
trapping scheme consists of a narrow channel of 2 jam 
width, and four adjacent half-circular electrostatic con- 
trol traps of diameter 10/im, located 1/xra away from the 
channel edges. The control traps are again loaded with a 
steady state Xd density (by, e.g., CW excitation) which 
can be tuned to control the energy profile along the chan- 
nel, denoted by E(f). The traps' location thus defines a 
trapping region along the cannel (marked by A), with a 
length of « 30/im. 
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FIG. 3. (color online). Upper panel: a proposed trapping 
scheme that consists of a narrow channel along the y direction, 
and 4 adjacent control traps. The effective trapping region, 
A, is marked by the shaded rectangle. Lower panels: (a) 
The energy profile along the channel with trapping on (green 
line), and off (blue line), (b) Initial Xd distribution in the 
channel at t = 0. The dashed lines mark the channel edges. 
The resulting Xd distribution for tf = 2/xs is shown (color 
scale x50) for (c) trapping "off' (empty control traps) and 
(d) trapping "on" (occupied control traps), showing a strong 
suppression of the Xd spreading along the channel. 

To study the dynamics of the proposed device, we use 
a numerical integration of the exciton hydrodynamics 
equation [H, 



dn 

~dt 



V (I D 



Jd + Je. 



n 

T 







(4) 



Here, n and r denote the (time and spatially depen- 
dent) Xd density and lifetime, respectively. The terms 
appearing inside the brackets are the diffusion induced 
current Jd = —DVn, the dipolar repulsion induced cur- 
rent Jd — — n/iaVn and a current due to external forces 
Jext = njiVE^ where \i is the exciton mobility and a is a 
material and structure parameter (see ref. [24]). We also 
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assumed a uniform temperature of the excitons through- 
out the sample. The principle of remote interactions is 
demonstrated by trapping freely expanding Xd inside the 
one dimensional channel, as described in Fig. [3] The 
double-bump energy profile along the channel is created 
by loading the control traps with Xd (Fig. G2a)). An 
initial gaussian distribution spot rii(f) = rii exp(— 4r 2 ) of 
Xd is created inside the channel at t = (Fig. [3f b)), 
and the system evolves in time according to Eq. ((4]) 
to tf = 2fis. The resulting Xd density profiles (Fig. 
C3(c,d)) show that the Xd are much better confined along 
the channel when the control traps are loaded (trapping 
"on"). 
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FIG. 4. (color online). The normalized Xd density, rif/m, at 
tf — 2 lis for changing (a) bump height Eb, (b) initial density 
rii and (c) temperature T. The trapping efficiency te as a 
function of the same parameters as in (a-c) is shown in (d- 
f). The dashed lines in (e) and (f) marks the density and 
temperature, respectively, where J ex t ~ Jd + Jd- 

To characterize and quantify the performance of such a 
flow control device, we analyzed the same device configu- 
ration described in Fig. [3j but under various conditions. 
Figure HJ^a-c) shows the resulting X^'s density profile, 
n/(r), along the channel, normalized to the maximal ini- 
tial Xd density n^, and its dependence on rii and the 
temperature T. One can see that rif increases with the 
bump height E^, and decreases with the initial X^'s den- 
sity rii and temperature T [25]. An objective measure for 
quantifying the performance of the device is the effective 
trapping efficiency (te), defined as 

Jn f (r)d 2 f 

^ = - — — ^ • exp(t / /r), (5) 
J rii(r)d r 

A 

where the integral is over the trapping area A (as marked 
in Fig. [3j) and the factor exp(t//r) compensates for the 
recombination loss of particles due to their finite life- 
time r, thus making the above definition of te a time- 
independent one (excluding other loss mechanisms such 



as boundary effects). The results show that there is an 
increase in te when E^ is raised (Fig. H^d)), and a de- 
crease in te with increasing rii ((Fig. ffl^e)) and increasing 
temperature (Fig. ID^f)). 

We can qualitatively understand these results by 
means of comparing the different Xd currents. For an 
energy bump with a typical width W, the induced 
current is estimated by: J ext = nLiVE « riifiEb/W, 
while for an initial Xd density profile n^(r) with typi- 
cal width L, the dipolar repulsion current will be: Jd — 
—niidSJri « fiarii 2 /L. The characteristic diffusion cur- 
rent, using the Einstein relation for a low density, non- 
degenerate Bose gas [23], is given by: Jd = —DVn « 
jikBTrii/ L. One expects that the trapping efficiency 
should fall whenever the size of the external trapping 
current J ext becomes smaller than the internal expand- 
ing currents Ji nt = Jd + Jd, or in other words when 
\Jext/Jint\ < 1- As is shown in our results, the be- 
havior of the analyzed system follows this rule-of-thumb 
fairly well. For the given experimental conditions (with 
L/W ~ 3), the dashed line in FigHJe) marks the initial 
density where J ext « Jd + Jd, and the dashed line in 
Fig[4jf) marks the temperature where J ext ~ Jd + Jd- 
There is a substantial decrease in te upon increasing den- 
sity (Fig. H(e)) and increasing temperature (Fig. H^f)) 
above these lines. 

These simple and intuitive results can be implemented 
in the design of devices for trapping and flow control 
of excitons. For example, one can build a flow device 
by applying a voltage gradient along a channel [7], and 
switch the flow on and off by remote dipole interactions 
with Xd fluids in the control traps, in a device similar to 
the one depicted in Fig. EJc). 

In summary, we propose a method for new dipolar 
exciton flow devices, which utilizes the remote interac- 
tions between spatially separated dipolar exciton fluids. 
It allows for a direct and model-independent density cal- 
ibration measurement, as well as for excitonic fluid-fluid 
control and manipulation. We give simple and intuitive 
rule-of-thumb estimates for the efficiency of such devices. 
We also show good agreement between the full numerical 
model predictions and our simple estimates. Such a con- 
cept is quite general and can possibly be implemented in 
different types of dipolar fluids and other device schemes. 
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and T = 1.5K. 



Zimmermann, [25] For varying Eb, we used rii 

In the case of varying the initial densities, we used Eb = 
W. Hansen, 0.125meV and T = 4K. For varying temperatures, we 

Wegscheider, set ra — 10 9 cm~ 2 and Eb — 0.125raeV. 



